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ABSTRACT 

We develop a general formalism for analysing parameter information from non- 
Gaussian cosmic fields. The method can be adapted to include the nonlinear effects in 
galaxy redshift surveys, weak lensing surveys and cosmic velocity field surveys as part 
of parameter estimation. It can also be used as a test of non-Gaussianity of the Cosmic 
Microwave Background. Generalising Maximum Likelihood analysis to second-order, 
we calculate the nonlinear Fisher Information matrix and likelihood surfaces in param- 
eter space. To this order we find that the information content is always increased by 
including nonlinearity. Our methods are applied to a realistic model of a galaxy red- 
shift survey, including nonlinear evolution, galaxy bias, shot-noise and redshift-space 
distortions to second-order. We find that including nonlinearities allows all of the de- 
generacies between parameters to be lifted. Marginalised parameter uncertainties of a 
few percent will then be obtainable using forthcoming galaxy redshift surveys. 
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1 INTRODUCTION 

The estimation of cosmological parameters from galaxy red- 
shift surveys, Cosmic Microwave Background (CMB) sur- 
veys, weak lensing surveys and velocity field surveys has 
become a major industry in cosmology, due to the data ex- 
plosion the field is currently undergoing. To date the ma- 
jority of the analysis has considered the problem of param- 
eter estimation from the large-scale, linear regime. However 
much of the important information on parameters lies on 
smaller scales, where nonlinear gravitational clustering has 
distorted the initial field. In the case of galaxy redshift sur- 
veys additional nonlinearity may arise form the bias rela- 
tionship between galaxies and matter, while in the case of 
the CMB the original pattern is still evident, allowing one to 
test the temperature fluctuations for evidence of primordial 
non- Gaussianity. 

The industry standard for parameter estimation in 
the linear regime is the multivariate Gaussian Likelihood 
method. Within this framework the uncertainty on parame- 
ters can be estimated from the second derivatives of the like- 
lihood function, the Fisher Information matrix. Tegmark, 
Taylor & Heavens (1997) were the first to study the Gaussian 
Fisher Information matrix as a general tool for characteris- 
ing the information content of cosmological surveys, while 
Heavens & Taylor (1997) demonstrated its importance for 
the design of galaxy redshift surveys. Tegmark (1997) used 
the Gaussian Fisher matrix to investigate the accuracy of 
parameter estimation from the Sloan Digital Sky Survey in 
some detail, but did not calculate the effects of nonlinear 



clustering, redshift-space distortions, or nonlinear biasing. 
Jungman et al. (1996) applied a Fisher analysis to the CMB, 
while Hu & Tegmark (1999) used it to estimate parameter 
uncertainties in weak lensing surveys. 

While likelihood parameter estimation in the linear, 
Gaussian regime has been extensively studied, until recently 
little work has been done on the non-Gaussian generalisa- 
tion. Amendola (1996) and Rocha et al. (2000) have ex- 
tended the likelihood methods to include non-Gaussian con- 
tributions by expanding the likelihood function in a series of 
cumulants. Amendola considered the effects of degradation 
of parameter uncertainties due to the non-Gaussian broad- 
ening of the likelihood function, while Rocha et al. applied 
their likelihood to estimate residual non-Gaussianity in the 
CMB. 

In the analysis of galaxy redshift surveys, Matarrese 
Verde & Heavens (1997) have investigated a Gaussian like- 
lihood estimator for the bispectrum as a tool for extracting 
the mean density parameter and galaxy bias from galaxy 
redshift surveys in the nonlinear regime. This builds on the 
work of Fry (1994) who suggested that the bispectrum could 
be used to break the linear degeneracy between the bias pa- 
rameter and the density parameter, Q m - Verde et al. (1998) 
developed this further to include redshift-space distortions. 
Scoccimarro et al. (1998) and Scoccimarro, Couchman & 
Frieman (1999) similarly investigated the bispectrum, in- 
cluding the effects of redshift-space distortions and galaxy 
bias. Scoccimarro et al. (2000) have recently estimated the 
bias parameters from the combined IRAS-QDOT and 2Jy 
and 1.2Jy redshift surveys. 
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In this paper we develop a new method for exploring the 
estimation of cosmological parameters from scales going into 
the nonlinear regime. We expand the general non-Gaussian 
multivariate Probability Distribution Function (PDF) and 
derive the non-Gaussian Fisher matrix. Since the most rel- 
evant non-Gaussianity will be due to nonlinear evolution of 
the density field, we use second-order perturbation theory to 
find the dependance of the higher-order moments on cosmo- 
logical parameters and estimate the Fisher matrix to second 
order. We illustrate these methods using a realistic model 
galaxy redshift survey, but note that they can be similarly 
applied to a 2-d weak lensing, or 3-d cosmic velocity survey. 

The paper is organised as follows. In Section 2 we de- 
scribe the expansion of the multivariate probability distri- 
bution function in both the discrete and continuum limits. 
In Section 3 we develop the formalism for calculating the 
Fisher information matrix in the non- Gaussian regime. We 
introduce a new function, the parameter entropy, to study 
the likelihood surface in parameter space in Section 4 and 
calculate a general expression for non- Gaussian fields. These 
methods are illustrated in Section 5 where we apply them 
to a model galaxy redshift survey, including the effects of 
shot-noise, redshift-space distortions and galaxy bias. We 
show how nonlinear effects change the estimation of param- 
eters in two surveys; the 2-degree Field galaxy redshift sur- 
vey (2dF; Colless 1996) and the Sloan Digital Sky Survey 
(SDSS; Gunn 1995). Finally, we present our conclusions in 
Section 5. Three appendices deal with technical issues. 



2 EXPANSION OF THE MULTIVARIATE PDF 
2.1 The discrete distribution 

Cosmological data sets can be modeled by a random field, 
4>{x), that is completely defined by its multivariate Proba- 
bility Distribution Function (PDF), p[<f>(x)]. The field <j>(x) 
can be discretised and treated like an iV-dimensional data 
vector, </>, whose components, (pi, may represent, for exam- 
ple, an array of pixel values or the amplitude of a set of 
harmonic modes. The PDF can be written as the Fourier 
transform of a characteristic function; 



p[<j>] = D[J]G[J]e 



J 



(1) 



where J D[J] = JY[ d N J/ {2tt) n is a multidimensional inte- 
gral. The characteristic function can be expressed by a set 
of cumulants, defined by the series 



lnG[J] = ■^ji<t>n-<l>in)cJj 1 :.Jj., 



(2) 



where there is an implicit summation over the j n , running 



from 1 to N. The 



,) c are the n order cumulants, 



or connected moments, of the field. 

If only the second cumulant exists the distribution is a 
Gaussian. If the initial distribution is Gaussian we can sep- 
arate out the n = 2 term and then deal with the rest of the 
terms in the series. Expanding out this second exponential 
to first order, keeping only terms to n — 3 and substituting 
back into equation (|l|) we find 



Pl<t>] 



x (l — ^ DijkJiJj Jk^j exp (iJm<t>m) 
JD[J] oq)[-ijiCyJj] 

exp {lJm<j>m) 



x 1 — —D 



3! J dfodfcdfa 



(3) 



where dj = 



bi4>j) c is the data covariance matrix and 
Dijk = {(fii&i&k) c is the data skewness. In the second line 
of equation (W|) we have used the substitution ij — d/d(j). 
Evaluating the Gaussian integral gives 



p[4>] = 



, In d 3 



exp [ - |frC^fo] 
v / 27rdetC 



(4) 



This procedure is a multivariate generalisation of the 
Edgeworth expansion for a single variable (Juszkiewicz et 
al. 1995, Bernardeau & Kofman 1995). Finally, performing 
the differentiation we find 



p\4>\ 



l 



-D 



6 

exp [• 



ijk 



C%a Cj b C kb <; 



V27rdetC 



(5) 



where Ar <:t ... in ] = l/n\(Ai 1 ...i n + perms) is the symmetriza- 
tion operator. This non-Gaussian multivariate distribution 
is the first major result of this paper. 

A well-known problem with truncating the moments, 
as we have done, is that the resulting distribution is not 
a well-defined probability distribution function. In fact the 
problem lies not in the truncation of moments, as the char- 
acteristic function is still well-defined, but in the approxi- 
mate calculation of the PDF itself in equation (^|). Hence 
the Edgeworth series, although properly normalised, can 
produce non-physical oscillations and negative values if the 
variance and skewness are pushed too far. This is a break- 
down of the series approximation, and one should take care 
not to force the PDF into this regime. 

Equation Q can be generalised to arbitrary order. The 
input for the nonlinear multivariate PDF are the higher or- 
der moments of the field. These can either be left arbitrary 
and measured from the data using the PDF as a likelihood 
function, or estimated by perturbation theory and used to 
estimate other parameters from the Bayesian relation 



£{0\x) = p(0)p(x\O). 



(6) 



where we have defined 6 = (Oi , 82, 9 m ) as the Tri- 
dimensional parameter vector. Usually we will assume a uni- 
form prior, p{6) = const. 

In developing the methodology it will be useful to ex- 
press the perturbed Gaussian distribution as 



p[cp] = p [cp](l - X) 



(7) 



where po[4>] is the underlying Gaussian distribution and 
X is a small perturbation. It is worth noting that some 
of the results of this paper are quite general and apply 
to all cases where the multivariate PDF is expressible in 
this way, for example the harmonic oscillator approach of 
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Rocha et al. (2000), or multivariate generalisations of other 
non- Gaussian fields. 



2.2 The continuum distribution 

In most calculations it is easier to consider the continuum 
distribution of the function <j>(x). This can be found by tak- 
ing the limit N —> oo and substituting — » V j d 3 k/(2n) 3 . 
From here on we choose to use a field embedded in 3 spa- 
tial dimensions, but the method is easily generalised to an 
arbitrary number of dimensions. In the continuum limit the 
underlying Gaussian distribution is 



pom = 



a 1 / 2 



exp 



-~ / d 3 x±d 3 x 2 < 



(Xi 



)C 1 (xi,X 2 )(j)(x 2 ) 



(8) 



where A is a normalisation factor. The expression for X, in 
the continuum case, containing the higher-order moments of 
the field 0(x), is given in Appendix A. The covariance and 
bivariance functions are, respectively 



C(x 1 ,X 2 ) = {()>(x 1 )(j>(x2)) c 
D(x 1 ,X 2 ,X 3 ) = (<f)(x 1 )(f){x2)4'{X3)) c 

and the inverse function C~ l {x\, x 2 ) is defined by 
J d 3 x C~ 1 (xi, x)C'(x,x 2 ) = 5d(xi — x 2 ). 



(9) 
(10) 

(11) 



3 THE NON-GAUSSIAN FISHER MATRIX 

To examine the information content available in nonlinear 
data it is useful to construct the Fisher information matrix 
for parameters. One issue of interest is the flow of informa- 
tion about parameter values during nonlinear evolution. In 
one respect nonlinear evolution destroys information about 
initial conditions. But since information has to be preserved, 
it must inevitably be transported up the hierarchy of cor- 
relation functions. However nonlinear growth itself depends 
on the detailed values of cosmological parameters. The ques- 
tion is does nonlinear evolution create or destroy parameter 
information? 

To explore these questions we can define three useful 
quantities: the entropy of the system, 



S^(C) 



dp In p, 



(12) 



where C = — In L is the average log-likelihood evaluated at 
the true parameter values; the average gradient of the log 
likelihood in parameter space, 



(d) = — I dp di In p = 0, 



(13) 



where di is the gradient in the i direction of parameter 
space, and is identically zero at the maximum likelihood 
values; and the Fisher Information matrix, the curvature of 
the likelihood surface about its maximum, which quantifies 
our knowledge of a set of parameters, 



dpdidjlnp= dp (diln p)(dj In p). (14) 



The parameter covariance matrix can be constructed from 
the Fisher matrix via 



T = (60*86) = -F~\ 



(15) 



where SO is the displacement from the maximum likelihood 
point, SO = 0-0'. 



3.1 Two-point correlations 

For the case of a Gaussian PDF we may derive a simple and 
familiar result for the Fisher information matrix. Since all 
higher order correlations are zero, X = 0. The Fisher matrix 
for Gaussian fields is therefore 



dpodidj In p . 



(16) 



The calculation is simplified if we choose the random field 
4>(x), to be the set of Fourier modes, S(k), where the covari- 
ance matrix is diagonal due to spatial invariance; 



C(*x,fca) = (5(ki)5(k 2 )) 

= (2Tv) 3 P(k)S D (k 1 + k 2 ). 

In this case the log-likelihood becomes 



2£o = 



d 3 k S(k)S*(k) 



V 



d A k 



lnP(fc) 



(17) 



(18) 



(2tt) 3 P(k) J (2tt) 3 

where the second term is the normalisation factor modulo 
an unimportant constant term. Here we have used the ma- 
trix relation In det C = Tr In C before transforming to the 
continuum limit. For a finite volume survey we have also 
used the approximation 5z>(0) = V, where V is the sur- 
vey volume. For a Gaussian distribution the Fisher matrix, 
equation (|l(]), becomes 

= ^ f-0^(9i^P(k))(d j lnP(k)). (19) 



A similar expression to equation ( |19| ) was previously ob- 
tained by Tegmark (1997), using a different derivation. 



3.2 Higher-order correlations 

The non- Gaussian regime may be explored via the perturbed 
multivariate PDF, equation (Q). The gradient of the log- 
likelihood, (&), is again identically zero, while the Fisher 
matrix becomes 



Tij = - J dp (1 - X)(d i hxpo)(dj]iipo) 

+2 J dpoXydfi lnpo - J dpoXiXj/il - X). (20) 

If the initial PDF is Gaussian, lnpo is an even function of 
variables and X, given by equation (B3), is odd. Expanding 
the last term this expression reduces to 



Tij = J-o,ij + (XiXj)c 



(21) 

■)g, denote 



which is accurate up to X 4 . The brackets, 
averaging over the Gaussian distribution. 

Again the analysis is simplified if we work in a Fourier 
representation, assuming spatial invariance. In this case the 
bivariance is 

D(k 1 ,k 2 ,k 3 ) = (S(ki)8(k 2 )S{k 3 ))c 

= (27r) 6 B(fci, fc 2 , k 3 )S D (k 1 +k 2 + fe 3 ), (22) 
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which selects only those combinations of wavevectors 
that form closed triangles. We will henceforth refer to 
B(k\,k 2 ,kz) as the bispectrum. 

With equations ( |l7| ) and (^2|) the perturbation term, 
X, reduces to 



X 



d% d% 
6 J (27r) 3 (27r) 3 



Q(k 1 ,k 2 )5(ki)S(k 2 )S(-ki - fc a ), (23) 



where 



<2(fci, k 2 



B{k 1 ,k 2 , -fci - fca) 
P(k 1 )P{k 2 )P(\k l + k 2 



(24) 



This term derives from the </> 3 -order term in the perturba- 
tion of the general PDF. The 0-order term vanishes under 
spatial invariance, since the bispectrum vanishes faster than 
the variance on large scales. This final constraint will al- 
ways be true for a realistic PDF since the Jensen inequalities 
(Kendall & Stuart, 1969) state that 



(x a ) < {y b ) a/t , 



(25) 



where x is an arbitrary random variables and a > b are real 
numbers. This implies that 



(26) 



<«5 2 } { ' ' 

so that the ratio of the bispectrum to the power will always 
vanish faster than the variance. 

The log-likelihood in the nonlinear regime is 

C = Co + X + ix 2 (27) 

where Co is the 2-point likelihood function, (p^), and X 
is given by (^). This expression can then be used to find 
the maximum likelihood values of parameters for systems in 
which the covariance matrix is close to diagonal. 

We calculate the {XiXj)G term in equation ( pl[ ) making 
use of the cumulant expansion theorem to expand the 6- 
point correlation function into its connected parts. For the 
case of a Gaussian average only the two point terms are 
non-zero so that 



{<5i<5 2 — 5g) = (8iS 2 ) c {d 3 S4,) ciS^Sa) c + perms. 



(28) 



Of these terms only six survive due to spatial invariance. 
Equation (pq) essentially states that to the order we are 
working to, the covariance of bispectrum triangles is zero. 
Matarrese et al. (1997), who use a Gaussian likelihood for 
the bispectrum, calculate the covariance matrix of the bis- 
pectrum to higher order and show that correlations between 
triangles appear through the pentaspectrum. 

The non-Gaussian Fisher matrix can subsequently be 
written, 



V f d 3 fci d 3 k 2 



+ 6 J (27r) 3 (27r)« 
xP(kx)P(ka)P(\ki + k 2 



t Qi(ki, fe 2 )Qj(fci, fe 2 ) 



(29) 



where the subscripts i and j denote differentiation with re- 
spect to the parameters. This is the second major result of 
this paper, and is independent of the nature of the non- 
Gaussianity, other than statistical spatial invariance of the 
field. 

With equations (||) and @ we can now address the 
question of information flow. Since the diagonal terms {X 2 )g 



are positive, the addition of the bispectrum term can only 
add information. This is because the dominant effect, to this 
order, is the additional information brought by the higher or- 
der moments. Nonlinear evolution will also change the shape 
of the PDF and degrade parameter information, but we see 
here that this is a higher-order effect. 



4 THE PARAMETER ENTROPY FUNCTION 

The Fisher Information matrix provides a compact way of 
determining both marginal and conditional errors on a given 
set of parameters. However, in some circumstances it is use- 
ful to have more information about the detailed shapes of 
likelihood surfaces in parameter space. This may be achieved 
by generalising the definition of the entropy from Section ^ 
so that 



8(9) = (C(9)) 



dp' In p, 



(30) 



where a prime denotes that a quantity is evaluated at its 
true (maximum likelihood) value. With this function one 
can map out the likelihood distribution of parameter space. 

An advantage of the parameter entropy is that it can be 
calculated directly from the covariance and bivariance func- 
tions, without having to take parameter derivatives. This 
may be a noisy process if the spectra are generated nu- 
merically, ie directly from a Boltzmann solver. Noise in the 
Fisher matrix can spuriously break parameter degeneracies, 
and the entropy is one way to avoid this. Noise may also 
prevent the Fisher matrix from being positive definite, and 
hence unphysical. Strong correlations between parameters 
can also make the Fisher matrix ill-conditioned and numer- 
ically difficult to invert. The entropy clearly contains more 
information about the parameters and the shape of likeli- 
hood contours, which may be important for investigating 
degeneracies among parameters. Finally, the parameter en- 
tropy does not assume the parameter likelihood distribution 
is Gaussian, which is implicit in the Fisher matrix. A dis- 
advantage of the entropy is that evaluating the marginal 
parameter uncertainties can be difficult since the full distri- 
bution of S must be mapped out around its maximum in 
parameter space. 

In the Gaussian case the entropy is given by 



So = 



dp' lnp - 



(31) 



In terms of the power spectrum this gives, up to an unim- 
portant constant, 



V fk 2 dk ( P'(fc) , 



(32) 



The nonlinear evolution of this quantity follows from the 
results of the last section. Expanding the non-Gaussian PDF 
as in equation (Q) gives 

S = S + J dp' X' lnpo — y dp (1 — X 1 ) ln(l — X) (33) 



Again if po is a Gaussian the second term in equation ( |33[ ) 
vanishes. The final log term can be expanded to second order 
in X, ln(l -X)&-X- X 2 /2, yielding 
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S = S -(XX') g + -(X 2 }g, 



(34) 



This expression again simplifies using a Fourier represen- 
tation and, with the results from the previous section, we 
find 



We = 7T 



V f d% d 3 k 2 „ 2 



6 J (2tt) 3 (2tt) 3 
xP'(|fei+fe 2 |) 



Q 2 (k 1 ,k 2 )P'(k 1 )P'(k 2 



(35) 



and 

{XX')e = jj ^3 (2^3 Q(ki,k 2 )B'(k 1 ,k 2 , -fc 1 -fc 2 ).(36) 

Equations @, @ and (51) are the third major new result 
of this paper. 

The parameter entropy evaluated at the maximum like- 
lihood point, tells us about the order in the system and can 
be written 



1 2 
S = S - -(X } G , 

where 

_ V fk 2 dk 



(37) 



(38) 



again up to an unimportant constant. In the case of gravi- 
tational instability, the power spectrum will always increase 
with time, so the linear entropy, So, will also increase. As the 
field is Gaussian distributed, this corresponds to the max- 
imum amount of disorder possible - that is the Gaussian 
requires the least information to specify it. 

The effect of the non-Gaussian term in equation $57\ ) 
always acts to decrease the entropy of the system, relative 
to the Gaussian case. While this may seem counter to our 
usual expectation that systems evolve towards higher en- 
tropy states, the effect of non-Gaussianity here introduces 
order, or structure, into the system, reducing the entropy. 



begin to break down for the case where triangles of the bis- 
pectrum lead to the mixing of information from both large 
and small scales (see section ^|) . 

5.1 Gravitational instability 

For a galaxy redshift survey a major source of non- 
Gaussianity comes from nonlinear evolution due to gravi- 
tational instability. For Gaussian initial conditions gravita- 
tional evolution may be modelled via perturbation theory. 
The skewness of the mass density field in perturbation the- 
ory is given by (Goroff et al. 1986) 



B(ki,k 2 ,k 3 



where 



(27r) 3 [2J(fe 1 , k 2 )P(ki)P(k 2 ) + cyc(23, 13)] 
xfo(fci +k 2 + k 3 ), (39) 



J(fei, k 2 ) = 1 + -^12 



k 2 + k 2 
k ± k 2 



«(On)[Mi2-l], (40) 



and where [i\ 2 — k\.k 2 . The dependance of J(k\, k 2 ) on the 
density parameter is weak and comes in only through the 
function «(0m) = D 2 /D 2 ~ — 3/7Sl m 2//63 , the ratio of first 
and second order growth functions (Peebles 1980, Bouchet et 
al. 1992). For simplicity we use the Einstein-de-Sitter value 
of k — -3/7 (Bouchet et al. 1992). 

The CDM-type model power spectrum we shall use 
takes the form 



A 2 (fc) 
where 



:<5| f (fc///) 4 fc n - 1 T 2 (fc,/ l ,n m ), 



A 2 (fc) = k 3 P(k)/2n 2 



(41) 



(42) 



and T(k, h, f2 m ) is the transfer function as given by Bardeen 
et al. (1986). The fiducial model we choose has parame- 
ters (fa, h, n m ,bi,b 2 ) = (5 x 10~ 5 , 0.65, 0.3, 0.002 1,0.5). 
For simplicity we assume spatial flatness. Figure hi shows 
the spectrum used for subsequent calculations. 



5 APPLICATION TO GALAXY REDSHIFT 
SURVEYS 

To illustrate the effects of including higher order statistics 
in the likelihood function it is useful to calculate the Fisher 
matrix and entropy for a simple analytic model. Since much 
can be derived from the nonlinear regime of galaxy redshift 
surveys, we choose this as our model, with application to 
the 2dF and Sloan Digital Sky Survey. We include sufficient 
detail in the model so that it possesses most of the properties 
associated with galaxy redshift surveys. The random field of 
interest is therefore the biased, redshift space distorted mass 
density contrast 5 s = Sn s /n 3 . 

Again, for simplicity, we elect to work in Fourier space. 
On large scales it is both more efficient and more accurate 
to exploit the natural symmetry of the survey and work in 
spherical harmonics (Heavens & Taylor 1995, Tadros et al. 
1999, Taylor et al. 2000). While this is important for lin- 
ear analysis, it is not necessarily so for the nonlinear analy- 
sis. Nonlinear effects dominate on smaller scales where the 
Fourier representation can replace spherical harmonics to a 
high degree of accuracy. However, this simplification may 



5.2 Choosing bispectrum triangles 

While our approach so far has been to retain as much in- 
formation as possible, in practice it may not be possible to 
include all shapes of triangles in fc-space due to CPU restric- 
tions. It therefore makes sense to identify and study a subset 
of shapes that maximise the information content. We discuss 
this issue in more detail in Section ^. For simplicity we con- 
sider here only two distinct shapes: equilateral triangles and 
degenerate triangles. Degenerate triangles are constructed 
from two parallel wavevectors plus a third that is aligned 
anti-parallel with magnitude equal to the sum of the other 
two. A convenient subset of the degenerate category, and 
the configuration we shall consider in this analysis, has the 
magnitude of the two smaller wavevectors being equal. Con- 
straints are made by inserting an appropriate dimensionless 
operator into equations ( pg| ) and (^J). For equilateral trian- 
gles we use: 

Aeq = j 8 D {ki - k)5 D {k 2 - k)5 D (k 3 - k); (43) 

whereas for the subset of degenerate (DG) triangles one pos- 
sible constraint is 
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Figure 1. Fiducial matter power spectrum for a model galaxy 
redshift survey. The linear and nonlinear power in real-space 
(solid lines) are plotted alongside the angle averaged redshift 
power spectrum (light solid line) with small-scale damping in- 
cluded. The shot-noise per mode for both of the model surveys 
are also plotted (see table [ij 2dF (dot-dash) and SDSS (dotted). 



S D (ki - k)5 D {k 2 - fc/2)fo(fc - k/2). (44) 



For the triangles of interest the bispectrum can be written 
B(k) = yP 2 (lt) Equilateral (45) 

B(k) = 4P 2 (fc/2) - P(k/2)P{k) Degenerate. (46) 

We plot the bispectra in Figure |^, where we have defined 
the function 



A|(fe) 



k 3 y/B{k) 



2tt 2 



(47) 



which is useful to compare with A (k). 

For an ideal, finite survey, in the absence of shot-noise 
and redshift-space distortions, the nonlinear Fisher matrix 
for equilateral triangles can conveniently be reduced to 



3~ i i 



k 2 dk 
2^" 



1 + 



■A\k) 



x [diln A 2 (k)][dj In A 2 (k)]. 



(48) 



Similar expressions for the parameter entropy and for the 
case of degenerate triangles can be found in appendix B. 

The effect of including nonlinear terms is two- fold. 
Firstly we see that the bispectrum introduces a second, pos- 
itive term into the Fisher matrix, driving the parameter un- 
certainty down. In addition, the nonlinear correction to the 
likelihood allows one to push the wave-range to higher k, 
again driving down the uncertainty. 
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Figure 2. Fiducial bispectra for the redshift survey model for 
equilateral triangles (bottom panel) and degenerate triangles (top 
panel). The real-space bispectra (heavy solid lines) are shown 
alongside the rcdshift-space bispectra (light solid lines) with small 
scale damping. The P(k)/n noise term (dotted) and the 1/n 2 
noise term (dot-dashed) are also shown for a 2dF-type survey. 



5.3 Poisson sampling of the density field 

Assuming the measured galaxy population can be modelled 
as a random sampling of a smooth underlying field, the 
power- and bi-spectra transform according to 

P(k) -> P(k) + -, 
n 

S(fei,fe 2 ,fe 3 ) -> B(fci,fc 2 ,fc 3 ) 

+ -[P(k 1 )+P(k 2 ) + P(k 3 )] + \, (49) 

where n is the mean number density of sources. If the noise 
varies across the survey, as is the case for flux-limited galaxy 
redshift surveys, this can be included in the models by the 
substitution n = n(r) and V — ► J d 3 r in the Fisher and en- 
tropy estimates. For simplicity we work with the volume lim- 
ited case and consider two models characterised by the num- 
ber density, n, and survey volume, V, as given in Table 111 
Noise levels for the surveys are plotted in Figures [j] and |2^ 
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Survey 


Mean Number 


Effective 




Density (n) 


Volume (V) 


2dF 


3 X icr 3 


9 x 10 7 


SDSS 


10 x icr 3 


10 x 10 7 



Table 1. Parameters for two model redshift surveys considered 
in this paper. 



5.4 Bias and redshift space distortions 

The use of redshift coordinates as radial distance introduces 
a distortion in the observed density field; 



S s (k) = D(kan)(l + I3^ 2 )8{k), 



(50) 



where the linear distortion modulates the density by a factor 
(1 + f3 n 2 ) (Kaiser 1987) with (3 = Sl^ 6 /6i and bi the linear 
bias parameter. The small-scale pairwise velocity dispersion 
can be modeled by convolution of the density field with a 
random component along the line of sight (Peacock 1992, 
Peacock & Dodds 1994); 



D(kafi) 



+ fcV/i 2 /2 



(51) 



where fi — f.k. Galaxy bias can be modeled by a local ex- 
pansion of the density field (Fry and Gaztanaga 1993), 



ml 



(52) 



where the coefficients b m are the bias parameters. The power 
spectrum for a biased, redshifted, and noisy galaxy survey 
is therefore 



6 RESULTS FOR A MODEL GALAXY 
REDSHIFT SURVEY 

6.1 Overview of parameter estimation 

In this Section we apply the model outlined in Section ^ 
to estimate parameter uncertainties from two model galaxy 
redshift surveys, based on the 2dF and SDSS (see Table 1 
for survey parameters). The parameters of the model are 
(h, 5b, fim, bi, 62), with fiducial values given in Section 5.1 
The model enters the analysis via the redshifted galaxy 
power spectrum, equation (JH^), and the bispectrum for 
equilateral triangles, equation (|55j), and degenerate trian- 
gles, equation (p6|) . These are used to estimate the function 
Q(ki, &2, £3), given by equation (pi|), which is then used to 
calculate the nonlinear Fisher matrix, given by equations 
( p9| ) and (|l9|). The parameter entrop y fu nction is calculated 
from equations @, @, @ and @. 

This section is laid out as follows. In Section 3.2 we 
present results for the marginalised parameter uncertainties 
for the two surveys, as a function of the truncation wavenum- 
ber for the Fourier integrals. Our basic result is that cosmo- 
logical parameters may be estimated to accuracies of a few 
percent using redshift surveys independantly of other mea- 
surements. In Section 6.3 we show the correlation coefficients 



for all the parameters, again as a function of wavenumber, 
and illustrate that to second-order all of the degeneracies 
between parameters are lifted. We illustrate how parame- 
ter estimation is i mpr oved using a simplified two-parameter 
model in Section 



6.4 



We explore the relative information 
con tent of equilateral and degenerate triangles in Section 
We begin with an estimation of the parameter uncer- 
tainties from the full model. 



6.5 



P°(k) = biD 2 (kafi)(l + PiffP{k) + -. (53) 

Note that we only consider the case of a deterministic bias, 
but this formalism can be easily generalised to include extra 
parameters for a stochastic biasing scheme (Pen 1998, Dekel 
& Lahav 1999). 

For equilateral triangles the angles between the line of 
sight and the wavevectors of the bispectrum triangles are 



Mi = M> 
M2 = 
M3 = -Mi 



1 

2^' 



M2, 



(54) 



where fii = r.ki. The corresponding redshifted bispectrum 
for equilateral triangles is (Hivon et al. 1995, Heavens et al. 
1998, Verde et al. 1998, Scoccimarro et al. 1999) 



B° 3 (k) 



D2(kafi) 



2b 3 1 Kei(k,n,(3,bi 



P 2 {k) 



+362(1 + /3 + 3/16/T) 
+-{D(km)(l + (3iA) 2 + cyc)P(fc) + \ (55) 

where i>2 is the second-order bias term and where kernel 
Ker(fc, /3, bi) is provided in Appendix B. Similar expres- 
sions for degenerate triangles can be obtained and are also 
given in the appendix. The redshifted power and bi-spectra 
are shown in Figures ^ and ^. 



6.2 Nonlinear parameter estimation from redshift 
surveys 

Marginal errors may be obtained from inversion of the Fisher 
matrix so that the the Icr error on parameter 0i is approxi- 
mately given by 



A9i = y/^' 1 ) 



(56) 



under the assumption that the likelihood distribution of pa- 
rameters is roughly Gaussian. In Figure ^ we show the re- 
sults of the Fisher analysis for the best attainable marginal 
uncertainties on the 5 independent parameters Sh, h, Q m ,bi 
and b2 as a function of k. Here, and in subsequent figures, 
k is the maximum mode, fc max , in all integrations. A 5-D 
joint estimation of parameter uncertainties with the Fisher 
matrix is only possible for a nonlinear likelihood since to lin- 
ear order the Fisher matrix is singular due to degeneracies 
between parameters. 

We give the results for two different redshift surveys 
which we define by the volume of the survey, V, and the 
mean number density of sources, n. Formally this represents 
a volume-limited survey, but here we shall use effective val- 
ues for a flux-limited sample. The volume and number den- 
sity are chosen to resemble those of the 2dF redshift survey 
and the Sloan Digital Sky Survey. The relevant parameters 
are found in Table ^. 

Interestingly, our results show that the accuracy of the 
parameter estimation is similar for the 2dF and the SDSS, 
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Figure 3. Marginalised error estimates on 5 parameters obtained from a Fisher analysis using the combined linear and nonlinear 
likelihood with equilateral and degenerate triangles. The results are shown for two sets of survey parameters resembling those of the 2dF 
(solid) and the SDSS (dot— dashed). Table 1 contains the survey parameters. 



for all of the parameters out to quite high k. The reason for 
this is that the effective volumes of the two surveys are sim- 
ilar (see Table 1), and the error scales like 1/VV. Although 
the SDSS has a factor of w 3 higher density of galaxies than 
the 2dF, this does not become important until the mod- 
els become noise dominated per mode. From Figures 1 and 
2, we see that for both surveys, this is not until around 
k — 1/iMpc -1 , well below the scale where our parameter 
estimation is valid. 

Figure ^ (LHS) shows the marginalised uncertainties on 
the parameters 5h, h and ft m . At k = 0.1 /iMpc -1 , normally 
the limit of linear analysis (Heavens & Taylor 1995, Tadros 
et al. 1999, Hamilton, Tegmark & Padmanabhan 2000), the 
predicted uncertainties lie around the « 30% level for both 
h and fl m . Estimates for 5h, the horizon scale amplitude 
of matter perturbations, fare a little better at around 20%. 
Extending the analysis to k — 0.3 feMpc -1 , near the limit of 
second-order perturbation theory, we find the uncertainty 
falls by an order of magnitude to « 2% for all three param- 
eters. 

Figure | (RHS) also shows the marginalised parameter 
uncertainty for the bias parameters bi and bi. For both sur- 
veys, truncated at k = 0.1/iMpc -1 , we find that bi is not 
detected at all, with bi detected at the as 10% level. How- 
ever, if the analysis is pushed to k = 0.3/iMpc -1 the error 
on bi drops to around a percent, while b2 is now measured to 
around 5% accuracy. While this looks encouraging for 2dF 
and the SDSS, we should caution that we have assumed a 
volume limited survey, with effective volume and noise levels 
chosen to mimic those of each survey. A more accurate as- 
sessment of the improved accuracy will require more detailed 
modeling and testing with N-body simulations, although the 



present study should provide a good indication of the mea- 
surement uncertainties. 



6.3 



Lifting all the parameter degeneracies 
redshift surveys 



in 



As well as the marginal uncertainties it is interesting to con- 
sider the correlations between parameters. The information 
about such correlations lies in the off-diagonal components 
of the Fisher matrix. Defining the correlation coefficient for 
parameters 9i and 9j as 



(57) 



we can quantify the degree of correlation as a function of 
maximum wavenumber. The results are plotted in Figure 



s 



for all 10 of the correlation coefficients, r^, in the 5- 
dimensional parameter space considered above. The indices 
on the correlation coefficients are given in Table || 

At low-fe the Fisher matrix is dominated by the linear 
contribution to the likelihood and subsequently all of the pa- 
rameters remain perfectly correlated/anti-correlated (ry = 
±1). The reason for this is that in the linear regime, one can 
only measure the three composite parameters, /3 = f2^; 6 /bi, 
from redshift-space distortions, biSn from the amplitude of 
linear galaxy clustering, and the spectral shape parameter, 
T = fl m h, which determines the break scale in the matter 
power spectrum. 

When the nonlinear terms begin to become important, 
at k — 0.04 /iMpc -1 , the Tij deviate from ±1 indicating 
the decoupling of the parameters. The reason for this uni- 
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Figure 4. Correlation coefficients for every pair of parameters in 
the 5-D marginalisation. The coefficients were calculated for the 
case of a 2dF-type survey and for a combined linear and nonlinear 
likelihood and both types of triangles. The key references Table ^ 
for the index of each parameter. 



versal removal of degeneracies lies in the lifting of the de- 
generacy between the bias parameters, b 1 and &2- With a 
single triangle shape, one can only measure the combina- 
tion 6i/(l + 62/2 J(fei,fe 2 )6i) (Matarrese et al 1997). With a 
second triangle shape this constraint is lifted. Figure [5] illus- 
trates this, where we have plotted the two degenerate like- 
lihood contours for bi and 62 for equilateral and degenerate 
triangles. The angle of intersection of these contours depends 
on the choice of triangles and is maximised by choosing very 
different shapes, motivating the choice of equilateral and de- 
generate triangles. With the bias degeneracy removed, the 
degeneracies in /3, F and &i& are all lifted. 

While the parameters are no longer degenerate, many 
remain quite tightly correlated. The amplitude 8h appears 
to become least correlated, becoming almost completely de- 
coupled from the other 4 parameters. This is because to 
linear order it is only degenerate with 61. 

As k gets larger, the rij coefficients begin to show de- 
viations from a steady decline; in some cases this results in 
them actually growing larger. This occurs because the cor- 
relations between parameters ultimately depend upon shape 
and size of the l-cr likelihood surface around its maximum. 
As this surface shrinks with increasing parameter informa- 
tion, it can change its shape, depending on the relative re- 
sponse to clustering information. This is what ultimately 
determines the correlation coefficient values. 



Figure 5. Lifting the degeneracy of the bias parameters using a 
combination of different triangle shapes for the bispectrum. The 
contours denote the 1, 2 and 3<r confidence regions for Gaussian 
entropy plus equilateral triangles and Gaussian plus degenerate 
triangles (faint lines) and for the combined effect (dark lines). 
The entropy has been marginalised over the amplitude <5jj. Red- 
shift distortions have been ignored and the noise and volume are 
similar to those of the 2dF 



index 


parameter 


1 




2 


h 


3 




4 


61 


5 


l>2 



Table 2. Cosmological parameters and their index for the corre- 
lation coefficient. 



6.4 A two-parameter model 

While the Fisher matrix results show the value of includ- 
ing nonlinear effects in parameter estimation, it is useful 
to examine in some more detail how nonlinearity increases 
parameter information. Rather than deal with all of the pa- 
rameters, we shall restrict our analysis to two parameters. 
One is the linear bias factor, 61, who's measurement allows 
us to lift all the parameter degeneracies. The second is the 
parameter combination, T = fl m h. The advantage of re- 
stricting our analysis to these two parameters is that their 
Fisher matrix in the linear regime is not singular, and can 
be compared directly with the nonlinear regime. We shall 
assume all other parameters are know. 

In Figure ^| (top panel) we have plotted the deriva- 
tives of the nonlinear, angle-averaged, redshifted function 
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Figure 6. Parameter derivatives (top panel) and 1<t error levels 
for the joint estimation of b\ and Y in redshift space (bottom 
panel). The top panel shows the angle averaged derivative of the 
Q function from equation ( p^ ) for equilateral triangles (solid lines) 
and degenerate triangles (dotted lines) taken with respect to each 
of the parameters. The bottom panel shows the joint errors on the 
parameters when using the Gaussian likelihood (solid line), when 
adding in the contribution from equilateral triangles (dashed line) 
and when adding in both equilateral and degenerate triangles 
(dotted line). The survey parameters closely resemble those of 
the 2dF, as listed in Table \ 



Q(ki, k2), given by equation (^) with respect to each of the 
parameters. This effectively shows the information content 
of the modes, with respect to each parameter. Dark lines 
show the derivative with respect to T, while lighter ones 
with respect to b 1 . Solid lines are for equilateral triangles, 
while dotted lines show the contribution from degenerate 
triangles. Figure ^ shows that degenerate triangles have a 
higher derivative, indicating that they contain more infor- 
mation that equilateral ones. We shall return to this point 





0.195 

r 



0.195 

r 




Figure 7. Entropy contours in the 6i — Y plane. The contour lev- 
els correspond approximately to 1 and 2<r confidence limits. From 
top left to bottom right the plots show: the Gaussian entropy; the 
Gaussian entropy (light lines) and the contribution from equi- 
lateral triangles (dark lines); the Gaussian entropy (light lines) 
and the contribution from degenerate triangles; and the Gaus- 
sian entropy (light lines) with the sum of the contribution from 
the Gaussian part plus equilateral and degenerate triangles (dark 
lines). For all plots k m ax was 0.2 hMpc - . The survey parameters 
closely resemble those of the 2dF. 



in Section 3.5. The similarity of the shapes of the derivative 
also suggests that these parameters are strongly correlated. 

In Figure |^ (bottom panel) we show the marginal la 
uncertainties on 6i and T obtained from the Fisher analy- 
sis. It is immediately striking that inclusion of just a sin- 
gle triangle configuration has a significant effect on the er- 
rors. Even integrating out to relatively modest wavenum- 
bers, k ~ 0.2 ZiMpc -1 , there may be a 20 - 30 % reduction 
in the marginal uncertainty. This effect continues to grow as 
more modes are included. Using both triangle configurations 
has a much greater effect, reducing the errors by nearly an 
order of magnitude out at k = 0.3 /iMpc" 1 . 

In Figure [j] we plot the entropy contours for a slice 
in the &i — Y plane, assuming the other parameters are 
known. The top-left plot shows the Gaussian contribution 
to the parameter entropy, where the maximum wavenumber 
is fc max =0.2 faMpc" 1 . The two parameters are highly corre- 
lated, and have uncertainties of A&i w 0.05 and Ar ~ 0.005. 
The top-right plot show the effect of adding nonlinear infor- 
mation from equilateral triangles; a « 30% improvement. 
By adding only degenerate triangles to the Gaussian term, 
in the bottom left figure, the uncertainty reduces more dra- 
matically. Finally, in the bottom-right plot we show the com- 
bined effect of all of the terms. 

Interestingly although the marginal errors are reduced, 
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inclusion of the nonlinear terms does not, in this case, make 
the two parameters, 61 and F, significantly less correlated, as 
suggested by the parameter derivatives in Figure Id. This can 
also be seen from the entropy contours in Figure m which are 
all rather elongated. The origin of this degeneracy is the ef- 
fect both parameters have on the power spectrum. The linear 
bias factor, 61, affects the amplitude of the power spectrum 
(ignoring its effects on redshift space distortions), while the 
shape parameter V shifts the break scale in wavenumber. If 
the power spectrum were a power-law, these effects would 
be degenerate. But since there is a break, this degeneracy 
is broken on scales below the break scale. Adding nonlin- 
ear information, while reducing the absolute error, does not 
reduce this correlation. 

6.5 The information content of triangles 

It is apparent from Figures ^| and ^| that degenerate triangles 
seem to contain a great deal more information about cos- 
mological parameters than equilaterals. As this result seems 
somewhat unintuitive, it is perhaps worth considering how 
the bispectrum triangles add information to the Fisher ma- 
trix. 

The effect of different triangles can be partly under- 
stood from consideration of equation For the non- 
Gaussian term, the Fisher matrix depends upon the deriva- 
tives of the function Q(k\,k,2) weighted by the product of 
P(ki)P(k2)P(\ki + fc2|). For equilateral triangles this fac- 
tor is the cube of the power. However, for triangles where 
the magnitudes of the wavevectors are not equal there is a 
mixing of power from different scales. This implies is that 
information contained in a given triangle configuration is 
related to the shape of the power spectrum. Although the 
parameter derivatives of Q(ki, £2) are larger for degenerate 
triangles, indicating a larger information content (see Fig- 
ure the power- weighting will ultimately determine which 
contributes the most information. 

For monotonically rising P(k), equilateral triangles are 
most significant because small scales (large k) contain most 
of the power. Conversely for power spectra that are decreas- 
ing functions of k, the maximum signal comes from large 
scales, hence for a given k max , triangles that have a small 
side will give the most weight and provide the most infor- 
mation. In the CDM case the greatest power lies around 
the break scale, fc*. The result is therefore an admixture 
of the two competing effects; longwards of k* the spectrum 
rises and equilateral triangles contain the most information, 
shortwards and it is triangles that contain a short side that 
win out. In the limit this means that for CDM there is a 
"maximal" triangle for parameter estimatio n th at is degen- 
erate, in the manner discussed in Section ( |5.2| ), with sides 
fc* , k m ax , kmax — k*. In practise, estimating parameters from 
degenerate triangles may be more difficult, as there will be 
fewer independent low-fc modes within finite survey volume. 



7 CONCLUSIONS 

We have presented a method for including higher-order mo- 
ments in the likelihood functions of cosmological fields in 
such a way that the parameter dependencies of the non- 
Gaussian terms may be used for estimating cosmological 



parameters. This non-Gaussian correction generalises like- 
lihood analysis for application to fields that either contain 
intrinsic non-Gaussianity, or have become non-Gaussian due 
to nonlinear gravitational evolution. Such a generalisation is 
fundamentally important for cosmology since most of the in- 
teresting fields show non-Gaussian properties. Our method 
has the advantage that we attempt to deal with the proba- 
bility distribution of the field explicitly, allowing the natu- 
ral combination of linear and nonlinear regimes. In this re- 
spect our work differs significantly from existing techniques 
that separate out the nonlinear regime and apply a Gaussian 
PDF to the higher-order moments themselves (Matarrcsc et 
al. 1997, Verde et al. 1998, Scoccimarro et al. 1999, Scocci- 
marro et al. 2000). 

We have presented a general formalism for calculating 
the non-Gaussian Fisher matrix and the shape of the likeli- 
hood around its maximum, the parameter entropy function. 
A central result of this paper has been to show that to lowest 
order, the dominant effect of adding nonlinearity, or more 
generally a non-Gaussianity, is to increase the the parameter 
information. This effect dominates over the degradation of 
parameter information due to the nonlinear evolution of the 
shape of the likelihood function. 

Applying our analysis of nonlinearity to a simple model 
for a galaxy redshift survey, including the effects of shot- 
noise, redshift-space distortions and galaxy bias, we have 
found that nonlinear effects may be extremely useful for 
placing tight constraints upon cosmological parameters. Of 
crucial importance is the fact that even at second-order, 
degeneracies can be broken so that all of the cosmolog- 
ically interesting parameters may be estimated indepen- 
dently. While going to higher-order may place tighter con- 
straints on the nonlinear bias function through estimation 
of the series coefficients b„, essentially redshift surveys of- 
fer up a great deal of their information without the need 
to turn to progressively more intricate perturbative calcula- 
tions. In addition, this analysis suggests that galaxy redshift 
surveys can be used in isolation to determine cosmological 
parameters. This greatly enhances their power as parameter 
estimation can be compared from independent surveys, such 
as the CMB, cosmic velocity fields and weak lensing surveys 
and combined further to reduce uncertainties. 

Our analysis relies upon a multivariate generalisation of 
the Edgeworth approximation for the probability distribu- 
tion function. As such one should be careful to determine the 
range of applicability of second-order perturbation theory, 
and the series expansion of the likelihood function. It is well 
known from studies of 1-point distributions (Bernardeau & 
Kofman 1995, Juszkiewicz et al. 1995 Gaztanaga et al. 1999, 
Taylor & Watts 2000, Watts & Taylor 2000) that the Edge- 
worth PDF has limitations - displaying unphysical features 
when the variance or skewness becomes too high. However, 
these effects do not automatically render the Edgeworth a 
bad approximation. Comparison with N-body simulations in 
redshift-space (Watts & Taylor 2000) has shown that peak 
of the distribution remains a good fit even when the ex- 
tremities begin to behave badly. Fortunately we have found 
that a good deal of information comes from the nonlinear 
terms at relatively low-fe and have shown that the degenera- 
cies between parameters are broken even for modest values 
of wavenumber. This is good news from a computational 
point of view as well since the number of modes available 



© 0000 RAS, MNRAS 000, 000-000 



12 A.N. Taylor & P.I.R. Watts 



for analysis grows like k . While this is beneficial in some 
respects, it can also lead to computational problems if the 
number of modes gets too high. In his case data compression 
methods may need to be incorporated (Taylor et al. 2000). 
In a forthcoming paper (Watts & Taylor, in preparation) 
we rigorously test our model on N-body simulations in or- 
der to accurately determine the limitations, and address the 
specific issues associated with application to real data sets. 



ACKF OWLEDGEMENTS 



PIRW thanks the PPARC for a research studentship. ANT is 
a PPARC Advanced Fellow. The authors thank Alan Heav- 
ens for useful discussion about the bispectrum. 



REFERENCES 

Amendola L., 1996, MNRAS, 283, 983 
Bardeen J.M., Bond J.R., Kaiser N., Szalay A.S., 1986, ApJ, 
304, 15 

Bernardeau F., Kofman L., 1995, ApJ, 443, 479 
Bond J.R., Efstathiou C, Tegmark M., 1997, MNRAS, 291, 
L33 

Bouchet F.R., Juszkeiwicz R., Colombi S., Pellat., 1992, 

ApJ, 394, L5 
Coles P., Jones B., 1991, MNRAS, 248, 1 
Colless M., 1998, In proc. 14th IAP meeting, Paris 
Dekel A., Lahav O., 1999, ApJ, 520, 24 
Folkes S., Ronen S., Price L, Lahav O., Colless M., Mad- 
dox S., Deeley K., Glazebrook K., Bland-Hawthorn J., 
Cannon R., Cole S., Collins C, Couch W., Driver S.P, 
Dalton C, Efstathiou C, Ellis R.S., Frenk C.S., Kaiser 
N., Lewis I., Lumsden S., Peacock J. A., Peterson B.A., 
Sutherland W., Taylor K., 1999, MNRAS, 308, 459 
Fry J.N., 1994, Phys Rev Lett, 73, 215 
Fry J.N., Gaztanaga E., 1993, ApJ, 413, 447 
Gaztanaga E., Fosalba P., Elizalde E., 2000, ApJ, 539, 522 
Goroff M., Grinstein B., Rey S., Wise M.B., 1986, ApJ, 311, 
6 

Gunn J.E., 1995, In proc. American Astronomical Society 
Meeting, 186, 4405 

Padmanabhan N., 2000, MN- 



1995, MNRAS, 275, 483 
1997, MNRAS, 290, 456 
Verde L., 1998, MNRAS, 



301, 



Juszkiewicz R., 1995, 



Hamilton A.J, Tegmark M. 

RAS, 317, L23 
Heavens A.F., Taylor A.N. 
Heavens A.F., Taylor A.N. 
Heavens A.F., Matarrese S. 

797 

Hivon E., Bouchet F.R., Colombi S. 

ApJ, 298, 643 
Hu W., Tegmark M., 1999, ApJ, 514, L65 
Jungman G, Kamionkowski M., Kosowsky A., Spergel D.N., 

1996, Phys Rev D, 54, 1332 
Juszkeiwicz R., Weinberg D.H., Amsterdamski P., Chodor- 

owski M., Bouchet F., 1995, ApJ, 442, 39 
Kaiser N., 1987, MNRAS, 227, 1 

Kendal M.G, Stuart A., 1969, "The Advanced Theory of 

Statistics, Vol. 2" , Griffin, London 
Ma C, Fry J.N., 2000, |astro-ph/0003345 
Maddox, S., 1998, In proc. ASP Conf. Ser. 146: The Young 



Universe: Galaxy Formation and Evolution at Interme- 
diate and High Redshift 
Matarrese S., Verde L., Heavens A.F., 1997, MNRAS, 290, 
651 

Peacock J.A., 1992, MNRAS, 258, 581 

Peacock J.A., Dodds S.J., 1994, MNRAS, 267, 1020 



Peacock J. A., Smith R.E., 2000, |astro-ph/0005010 
Peebles P. J., 1980, "Large-Scale Structure in the Universe", 

Princeton University Press, Princeton 
Pen U., 1998, ApJ, 504, 651 



Rocha G, Magueijo J., Hobson M., Lasenby A., 2000, astro- 



ph/0008070 

Scoccimarro R., Colombi S., Fry J.N., Frieman J. A., Hivon 

E., Melott A., 1998, ApJ, 496, 586 
Scoccimarro R., Couchman H.M.P., Frieman J., 1999, ApJ, 

517, 531 

Scoc cimarro R., Feldm an H.A., Fry J.N., Frieman J., 2000 
festro-ph/0004087l) 



Seljak U., 2000, |astro-ph/0001493| 



Tadros H., Ballinger W.E., Taylor A.N., Heavens A.F., Efs- 
tathiou G, Saunders W., Frenk C.S., Keeble O., McMa- 
hon R., Maddox S.J., Oliver S., Rowan-Robinson M., 
Sutherland W.J., White S.D.M., 1999, MNRAS, 305, 
527 

Taylor A.N., Watts P.I.R., 2000, MNRAS, 314, 92 

Tayl or A.N., Ballinger W.E., Heavens A. F., Tadros H., 2000, 

astro-p h/0007048| 
Tegmark M., 1997, Phys Rev Lett, 79, 3806 
Tegmark M., Taylor A.N., Heavens A.F., 1997, MNRAS, 

480, 22 

Verde L., Heavens A.F., Matarrese S., Moscardini L., 1998, 
MNRAS, 300, 747 



Watts P.I.R. , Taylor A.N., 2000, astro-ph/0006192 



© 0000 RAS, MNRAS 000, 000-000 



Nonlinear Parameter Information 13 



APPENDIX A 



For completeness, we provide the general continuum form for the non-Gaussian perturbation term appearing in equation 

4 



X = jY\_d 3 x. l D(x 1 ,x 2 ,x 3 )C 1 (x 4 ,x 1 )C 1 (x 2 , x 3 ) 4>{xa) + cyc(231, 312) 



+ 



i = i 

6 



\[ d 3 x t D( 



Xl, X2,X 



z)C 1 (X1,X4)C 1 (X2,X 5 )C 1 (x3,xe) (j)(x4)<j)(xs)<p(x6) 



(58) 



where eye means to permutate indices. 



APPENDIX B 



In the absence of noise and redshift distortions the nonlinear Fisher matrix and entropy can be written in a particularly 
pleasing form. The Fisher matrix for equilateral triangles reduces to 



_ v f Pdk r 



+ (t)^ 



[& In A 2 (&)][&/ In A 2 (fc)], 



(59) 



while the parameter entropy is given by 

V fk 2 dk f, A2/1X . A' 2 (fc) 1 /12\ 2 . r/ A ' 2 (fc)\ 2 2A' 2 (fc) 



2tt 2 



A 2 (k) 



For the case of degenerate triangles the Fisher matrix becomes 

8 (A 2 (fc/2)) 2 



V f k 2 dk 
2 2tv 2 ' 



3 A 2 (fc) 

-^A 2 (fc/2) (9 [8 lnA 2 (fc/2) ft] In A 2 (fc/2)) } 
and the parameter entropy is 



[ft In A 2 (k)] [ft In A 2 (k)] + \ A 2 (k) [ft In A 2 (Jfe/2)] [ft In A 2 (fc/2)] 



r f '' 1,lk f lnA 2 (fc) I A ' 2{k) I 1 ±' 2 (k/2)\ A ' 2{k/2)A ' 2{k) 1A ' 2 ( fc /2)A' 2 (fc) , 1 A' 2 (fc/2)A' 2 (fc) 



2tt 2 



A 2 (fc) 12 



[A 2 (fc)] 2 



+ 



2 A 2 (fc/2)A 2 (fc) 16 [A 2 (fc/2)] 



lA' 2 (fc/2) 1 A' 2 (fc) + lA' 2 (fc/2) i lA' 2 (k) 



2 A 2 (fc) 32A 2 (fc/2) 8 A 2 (fc/2) 8 A 2 (fc) 



(60) 



(61) 



(62) 



APPENDIX C: REDSHIFTED BISPECTRUM 



Here we present some of the formulae used to calculate the redshift space corrections for the bispectrum. The kernel from 
equation (jHHJ) for equilateral triangles is, 

'2 !„,„,_ , , . 3 o2 



Ker 



3 {- + -0(9/7 + h) + — /3 2 [3 - 76i(2 - 61)] + — [3 3 [ M 2 (3 - 4 M 2 ) 2 + 76i(l - b?)(9(l - M 2 ) + 8 M 4 (3 - 2 M 2 )) 

--Ls 4 (&i - 1)(27- lUfi 2 + 384 M 4 - 256m 6 )}. 
256 J 



(63) 



The nonlinear small scale damping factor is 
1 1 



D2(kan) — 



y/1 + k 2 a 2 14/2 y / l + k 2 a 2 Li 2 /2 y/l + k 2 a 2 fi 2 /2 
1 

■sjl + 3/4fc 2 o- 2 + 9/64fc 4 a 4 + fc 6 crV(3 - 4^ 2 ) 2 



(64) 



The calculation in Section 5.4 can be repeated for degenerate triangles. In this case the angles between the line of sight 
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(65) 



and the wavevectors of the bispectrum triangles are 

A*i = 

= -A* 

(J,3 = 

The redshifted bispectrum is then given by 

B s g (k) = D 2 (k,ix) [2b?Ker(fc, fi,/3, bi) + 6 2 6?(1 + /3 Al 2 ) 2 (2P(fc) + P(fc/2))1 P(fc/2) 

+ i6 2 J D 1 (fc, Atl )(l + l3p,\) 2 P(k) + -blD^k,^! + f3(i 2 2 j 2 P(k/2) + -blD^k, fj, 3 ){l + f3n%) 2 P{k/2) + \, (66) 

71, 72 72 ft 

where in this case the kernel is 



Kcr 



4P(fc/2) -P{k) 



4(3 + 6i)P(fc/2) -3P(fc) 



(46i - 46? - 3)P(fc) + 2(56i + 6? + 6)P(fe/2) 



7/ J 



(1 - 861 + 8b 2 )P(fc) - 4(1 + 6? )P(fc/2) 



/9 3 - 61 At 6 



(2-26 1 )P(fc) + (l + 6i)P(fc/2) 



and the damping term 

£>2(*,/i) = 



+ fc 2 (7 2 A t 2 /8)(l + ifcVVTS) 



(67) 



(68) 
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